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The problem of vortex pair motion in two-dimensional plane radial flow is solved. 



in 

(N 



Under certain conditions for flow parameters, the vortex pair can reverse its motion 
within a bounded region. The vortex-pair translational velocity decreases or in- 



creases after passing through the source/sink region, depending on whether the flow 
is diverging or converging, respectively. The rotational motion of two corotating vor- 
texes in a quiescent environment transforms into motion along a logarithmic spiral 

(J 



in the presence of radial flow. The problem may have applications in astrophysics 



and geophysics. 



> ' I. INTRODUCTION 

o 



In the number of astrophysical and geophysical problems the situation arises in which a 



system of vortices (in the simplest case, a vortex pair) interacts with diverging or converging 
radial flow. For example, in the dipolar vortex model of an active galactic nucleus (AGN) 
[1], the obscuring tori 1 [2 - 5] are toroidal vortices in which selfgravitation is balanced by 
^ ■ centrifugal forces [6]. (Possible geophysical applications are discussed at the end of this paper.) 

b : 

The vortex motion inside an obscuring torus transform it from a pure geometrical object 
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1 The idea of obscuring tori underlies the unified model of AGN [3, 4] in which the difference between types of 
active nuclei, such as radio galaxies and quasars (as well as between types of Seyfert galaxies), is attributed 
to the position of the obscuring tori relative to the line of sight. For example, in the case of quasars is 
possible to observe the central AGN region adjacent to the black hole at the center of a galaxy, whereas the 
central region of a radio galaxy is obscured by the surface of a torus. In 2004, the idea was substantiated 
by direct observations of the obscuring torus with an optical interferometer combining 8.2 m telescopes at 
the Southern European Observatory [5]. 
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FIG. 1: Dipolar vortex model of AGN obscuring tori (symmetry plane cross section) [1]: M c is the 
central mass and R is the major radius of a torus. 

to the dynamical one that allows to describe dynamical processes in AGNs. The vortex motion 
can be caused by interaction between the torus and the radial outflow of wind and radiation 
from the AGN central region, which is responsible also for the dipolar structure of the vortex 
motion in the torus (see Fig. 1). 

The foregoing discussion motivates a study of the dynamics of interaction between vortices 
and radial flow. Problems of this kind seem to have never been analyzed. The simplest model 
is the planar interaction between radial flow and a dipole vortex pair. (The corresponding 
flow pattern can be envisaged from Fig. 1, but the present analysis is restricted to a plane 
two-dimensional problem.) In this paper, we present an exact solution for the vortex motion in 
an inviscid fluid. Possible exactly solvable generalization for rotating radial flow is considered 
at the end of the paper. 



II. PLANAR VORTICES IN DIVERGING RADIAL FLOW 

The input time-dependent streamfunction equation in the plane case has the form [7] 

dAib 

^f + J(iP,AiP) = 0, (1) 
where the flow velocity v is related to the streamfunction ip as 



dip dip 
dy' v dx 



v x = -—, v y = — , Aip = curl z v 



and J(a, f3) is the Jacobian determinant. The equation (1) is the z component of the vorticity 
conservation equation in two-dimensional case 

d , „ d d „ 
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We represent the streamfunction as a decomposition into a regular part ip r describing the 
flow and a singular part ip s whose components describe point vortices: 



(3) 



where 



A m is the mth vortex strength, and r m is the corresponding position vector. 2 
The singular component satisfies the Poisson Equation [9, 10] 

A^ s = ^2 A mS(x - x m )5{y - y m ). 

m 

Substituting (3) into (1) we obtain the streamfunction equation 

dAip r 



(4) 



dt 



+ J(lp s + lp r , Mr) = 



(5) 



and expressions for the velocity components of the mth vortex 



Vr, 



dy 

d(A + ^) 



dx 



(6) 



where tp™ is the streamfunction ip s minus the contribution of the mth vortex. 

If the regular component ip r corresponds to a diverging radial flow with source at the 



origin, 

Vv — —QVi Q = const > 
then vorticity vanishes outside the source: 

Aip - | 1 d ^ r | 1 ^tl - Q 

dr 2 r dr r 2 d(f 2 

Now, Eq.(5) becomes an identity, and expressions (6) combined with (7) yield 



(7) 



(8) 



x r 



dipl 



Vm = 



dy 
OA 



x 



dx 



(9) 



2 We use the approach developed in [8] . 

3 Note that the assumption of incompressible flow is inapplicable within a small neighborhood of the radial 
flow source depending on sound velocity. 
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(10) 



In the case of two vortices (Fig. 2), Eqs. (9) reduce to 

A 2 yn , 

■'2 - — -•- (>— ■ m = 

where (xi2,yn) = {x\ — X2,Ui ~ IJ2) and 7*12 = |ri — r 2 |. Consider a dipole vortex pair 

yt 




■J 




FIG. 2: Two point vortices, 1 and 2, 
in Cartesian and polar coordinates. Ar- 
rows directions correspond to A\ < and 
A 2 > 0. 



FIG. 3: Vortex pair in diverging radial 
flow. 



(A2 = —Ai) symmetric about the x axis (Fig.3) in the initial moment of time. Since X12 = 
and y 2 = —yi, it will suffice to analyze the equations of motion of the one vortex located at 
y = yi > 0: 

A x 

x = -- — + Q— , (11a) 



Airy r 

y 



y = Q 



(lib) 



where x = x\ = x%, r = x 2 + y 2 , and A = A2 > 0. In a quiescent environment, the pair 
moves with velocity A/4iry antiparallel to the x axis. It follows from Eq.(llb) that y > i.e., 
the separation between the vortices increases with time elapsed, while the interaction between 
them (the first term in Eq.(lla)) decreases in a diverging flow. 
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III. HAMILTONIAN FORMULATION AND INTEGRATION OF GOVERNING 

EQUATIONS 

An arbitrary system of vortices satisfying Eqs.(9) can be described in terms of the Hamil- 
tonian 

H = — y~] A m A n In r mn - Q ^ A m arccot— . 

rn^n m 

Indeed, multiplying Eqs.(9) by A rn and verifying that 

. dH . dH 

Xm^ = 0, 771 = 1,2,..., 

dx m dy m 

we can write 

dH dH 
oy m ox m 
In the case of a symmetrically moving vortex pair, Eqs.(ll) can be rewritten in Cartesian 

coordinates in the canonical form 

dH . dH 

x = -w y= ^- (12) 

with the time-independent Hamiltonian 

A x 
H(x, y) = — lny — Q arccot— . (13) 

The trajectory of the vortex pair is parameterized by the "energy integral" of the system (11): 

A x 
E = — \ny - Q arccot- . (14) 
Ait " y 

Hence, the "phase" of movement along the trajectory $ = arccot (x/y) can be expressed as 

A , E 

^ = ^Q lny -Q- (15) 

Depending on vortex and flow parameters relation the case of reverse vortex motion may 
be possible. The vortex pair reverses its motion between x_ and x+, where 



■2 -Q / / / .1 ? 



x ± = ^|l±,/l-^j | (16) 

are the turning points where x changes sign. For existence of the turning points it has to be 
fulfilled the inequality 

(2nQ/A) 2 > 1 . (17) 
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If A > 0, then the pair moves antiparallel to the x axis and infinity corresponds the phase 
values 

$ = +0 , x = +oo , 

(18) 

$ = 7T — , X = —OO . 

If A < the pair moves parallel to the x axis. The cotangent argument lies in this case in the 
interval (—it, 0) and the phases corresponding to infinity are 

$ = -0 , x — -oo , 

(19) 

$ = -7T + , X = +00 . 

It convenient to regard the vortex movement in polar coordinates, x = rcosip, y = 
r sirup, ip = arccot(x/y). In this case the Eqs.(ll) take the form 

r = --^-coty? + — , (20a) 
4nr r 

These equations can also be represented in Hamiltonian form: 

dH . dH 

(21) 



A A 
H = — lnf + — In sin 2 ip - 2Qip , 

4lX 4lX 

where £ = r 2 and Hamiltonian H(£, (p) is twice as large than Hamiltonian H(x, y) in Cartesian 
coordinates H(£,<p) = 2H(x,y) for the same system. Thus, the Hamiltonian equations of 
motion written in these coordinates are 

A 

£ = 2Q - — cot V> 

27T _ A (22) 

The distance from the center of source to the vortex reaches a minimum r* when = 
arccot(47rQ/A). The equation of trajectory r(ip) parameterized by the energy integral is 
solved directly relative to the distance from the center: 



£(<p) = r\<p) = exp \^—E + -^y? " In sin 2 y?J (23) 

where E = 2E and E is given by (14). Using (23), we integrate Eq.(20b) in quadratures: 

4tt ft 

t — to — — I d(p i exp (2(/i + X(pi) — In sin 2 ip{j , (24) 
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where 



AtiE 



A = 



4vrQ 



A ' A 

The function t(tp) may be expressed in terms of elementary functions. The inversion of this 
function and calculation of r(t) can be performed numerically. 

Changing back to the Cartesian coordinates, we obtain the trajectory equation 



x 



y cot 



A E 



(25) 



Hence, the vortex locations at infinite distance are found by using (19): 

[A-kE s 



y- 



y+oo 

- exp 



= exp v A 
An(E + ttQ) 



(26) 



Thus, the vortex separation increases by a factor of exp(4n 2 Q/A) after the pair has passed 
through the flow region. 

IV. BLOW-OFF OF THE VORTEX PAIR COMPONENTS 



Suppose that the energy corresponds to the desired trajectory. For example, condition 
(17) must be satisfied for a trajectory with reverse motion. Using the expression for the 
trajectory, we represent (16) as 



2ixQ 



x± = y± 
where the ordinates 

y± = exp 

are expressed in terms of the phases 

2vrQ 



1 ± 




— (E + Q- arccot^ ± ) 



(27) 



(28) 



i± Wi 



A 
4^Q 



(29) 



We take a value yo in the interval defined by (26) that lies sufficiently close to a turning point. 
Using the expression for the trajectory, we determine the corresponding value Xq, numerically 
calculate the function x(y), and invert it to find the required y(x). 

Figure 4 demonstrates that a vortex pair components moves apart as it approaches the 
source from infinity and can indeed reverse its motion along the x axis. Having traveled to 
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a distance sufficiently far away from the source, where the radial flow is weaker, the pair 
components resume their motion in the direction determined by the sense of vortex rotation. 
The separation at negative infinity given by (26) is reached asymptotically. 




FIG. 4: Trajectory of a vortex in the region 
of reverse motion. The arrow indicates the 
direction of motion of the vortex. The pa- 
rameter values Q = 1, A = 2ir — 0.6, and 
E = correspond to the trajectory passing 
through the point with xq = 6.17 and yo = 2. 



FIG. 5: Vortex trajectories for several values 
of L = 2ttQ/A and E = 0. The curve cor- 
responding to L = 1.0 separates the regions 
of reverse vortex motion [L > 1) and motion 
without turning points (L < 1). 



Vortex trajectories in radial flow can also be calculated by solving Eqs. (11) numerically 
with the values xq and yo above used as initial conditions (see caption to Fig. 4). Figure 5 
shows examples of trajectories with and without reverse motion. 



V. ACCELERATION OF A VORTEX PAIR BY CONVERGING RADIAL FLOW 

Consider a vortex pair in a converging flow, as shown in Fig. 6. Since Q < 0, we change 
to a new parameter according to Q — —P = const < and rewrite the expression for the 
regular streamfunction component as 



tpr = -Q<P = P<£ ■ 



(30) 



It is clear from Eq. (lib) that y < 0; i.e., the vortex separation decreases with time elapsed. 
According to (18), the separation at infinite distance is 



V+c 



exp 



AttE 



exp 



4n(E - ttP) 
A 



(31) 
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FIG. 6: Vortex pair in converging flow FIG. 7: Trajectory of a vortex pair compo- 

(in Cartesian coordinates). nent in the region of reverse motion. The 

arrow indicates the direction of vortex mo- 
tion. The parameter values are P = —Q = 1, 
A = 2vr - 0.6, and E = 0. 

It follows from (31) that the vortex separation decreases by a factor of exp(47r 2 P/v4) after the 
pair has passed through the sink region. 

The trajectory of the pair is found by following an analogy with the analysis above. 
It is clear that the vortex components separation decreases, and the pair can reverse its 
motion within a certain portion of its trajectory near the sink if the parameters satisfy certain 
conditions, as illustrated by Fig. 7. Having approached the sink to a sufficiently close distance, 
the pair accelerates as the vortex separation decreases and continues to move with a higher 
velocity determined by the sense of vortex rotation after it has passed the second turning 
point. At negative infinity, the pair reaches an asymptotic steady state with a smaller vortex 
separation and a higher translational velocity. In the example illustrated by Fig. 7, the velocity 
of the pair "thrown" out of the sink region is higher by three orders of magnitude. 

VI. LAGRANGIAN DESCRIPTION 

Hamiltonian (21) in polar coordinates can be represented as 



H = K(0 + U(<p) . 



Here, (p and £ are interpreted as generalized coordinate and momentum, respectively (cf.(22)); 
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if (£) = (A/ An) ln£ is the kinetic energy and potential energy is expressed as 

A 

U (ip) — — In sin 2 <p — 2Qp . 
An 

The corresponding Lagrangian is 

A A 

L((p, p) = — ln<p — — In sin 2 p + 2Qp + const , (32) 
An An 

where the generalized velocity ip = A/An^ is inversely proportional to £. 4 Given L(ip,ip), we 
integrate the Euler-Lagrange equation 

ip_ _ An dU(if) 
ip 2 A dcp 

to find the solution 



(33) 



dp ( 8nQp 



t-t = C -^exp i—^- , (34) 
J<po sm ip \ A J 

which is equivalent to solution (24) obtained by using the Hamiltonian approach. In this case 

„ An (%nE 
C = — exp 



A \ A 

VII. COROTATING VORTICES IN RADIAL FLOW 

We consider two corotating vortices with equal circulations symmetric with respect to the 
source of the radial flow (Fig. 8). 

It is well known [9] that the separation 2r between the vortices in a quiescent environment 
remains constant and they rotate about the center of vorticity with an angular velocity of 
A/2nr 2 . In radial flow, equations of motion (11) for a vortex are modified as follows: 

* = -^ + G£> (35a) 

y = ~ + Q ]L 2 - (35b) 

An r z r z 

where y — yi — —y2 and A± = A 2 = A. 

System (35) yields r 2 = 2Qt + r 2 ; i.e., the squared separation between vortices in a 
diverging flow is a linear increasing function of time (Fig.9a). In a converging flow, the vortex 



4 The inverse proportion between generalized velocity ip and momentum £ implies that with accuracy to 
unessential constant summand L = £tp — H = —H. 
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FIG. 8: Initial state of two corotating vortices with A± = Ai 
vortices circumrotate. 



A. In a quiescent environment, the 
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FIG. 9: Trajectory of a component of the pair of corotating vortices of equal strength A = Anx: 
(a) diverging radial flow with Q = 0.1; (b) converging flow with P = —Q = 0.1. 

separation vanishes in the finite time t = Tq/2P (Fig. 9b), while the angular velocity of vortex 
rotation indefinitely increases (see Footnote 3). 

Changing to the complex coordinate w = x + iy, we have |u>| 2 = r 2 . Accordingly, Eqs. 
(35) are equivalent to 

w /_ A \ 1 / A A I 

w \ An J \w\ A 

Integrating Eq. (36), we obtain 



4tt J 2Qt + r, 



(36) 



2Qt\ K 1+i 4^y 

W = Wq [ H 2" 



(37) 
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where wq is the initial value of the complex coordinate. In the limit of Q — > 0, the solution 
describes uniform rotation. Indeed, changing to the new variable r\ = 2Qt/r^, we have 



w 2 = w l . [(1 + rf/vyAt/inr* _ (3g) 



Since lim(l + r]) 1 ^ = e, we obtain 

?7 — >0 



w — > Wq exp % 



At 
2^1 



In the polar coordinates, the motion is particularly simple. It is obvious that the instantaneous 
angular velocity is expressed as follows (H = A ln£/47r — 2Q(p): 

A 



<p(t) = 



4n(2Qt + r 2 ) ' 
Hence, 



4 m + ^ 

SttQ V ^ 

i.e., the trajectory is a logarithmic spiral. The limit as Q — * is easily taken. 

VIII. VORTICES IN ROTATING RADIAL FLOW 

Flow rotation can be taken into account by adding a corresponding term to the expression 
for the regular streamfunction component, 

i) r = -Qip + ip? , 

where ^? = Qr 2 /2 represents rigid-body rotation with angular velocity Q. Then, Eqs. (10) 
are rewritten in terms of w = x + iy as 

dw 2 iA x Q ldyj 

_-,n W2 = -^(«, 1 - W2 ) + - ? «, 2 . 

In the rotating reference frame we substitute w = w'exp(ifit) and use the relations r 2 = 
|u>| 2 = \w'\ 2 and r 2 2 = \wi — u>2| 2 = — ^2 1 2 to derive the equations of motion for vortices 
in irrotational radial flow, which makes it possible to use the results obtained above. 
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IX. CONCLUSIONS 



We have shown that the problem of vortex pair motion in two-dimensional radial flow has 
an exact solution. Under certain conditions for flow parameters, the vortex pair can reverse 
its motion. In a diverging flow, the separation between the vortices increases, but the pair 
does not break apart. In a converging flow, the vortices move inwards, and the translational 
velocity of the pair increases. In both cases, the vortex separation in the pair changes by a finite 
amount as it passes through the source region. The squared separation between corotating 
vortices of equal strength is a linear function of time. In a diverging flow, they move apart 
indefinitely. In a converging flow, the vortex separation vanishes in a finite time. Note that the 
problem of vortex motion in rotating radial flow may have qualitative applications in planetary 
geophysics. For example, the recently discovered double-lobed vortex swirling around Venus 
south pole [12] may be described by a similar model if there is meridional flow in the polar 
region. The present model cannot describe vortices spun up by the wind, since circulation 
is conserved. Varying circulation can be treated by allowing for dissipation. It would be 
interesting to analyze the asymmetrical motion of a vortex pair in a radial flow. Finally, also 
of interest is the motion of vortex rings in radial flow. Translated by A. Betev 
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